x=0:1:20;y=0:1:15
[X Y]=meshgrid(x,y);
Q=atan((1609621605679839*Y)/4611686018427387904);
T =(1609621605679839*X)/31709018399331233792;
U=(884898835882229292997804032*atan((1609621605679839*Y)/4611686018427387904))./((52056435699881628425312815446547*Y)+eps);
 
surf(X,Y,U)